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I. INTRODUCTION 



The spectacular observations of Bosc-Einstein condensation (BEC) in atomic vapors reported last year Jlj^ have 
opened up new avenues of research into the physical properties and nature of Bose-condensed systems. Recent 
experimental developments include reports of a new trap capable of holding larger number of atoms and measurements 
of the condensate fraction and mean-field energy || , non-direct observation of the development of the condensate [H , 
measurements of the collective oscillations of the condensate ||- jjj, and an output coupler for an atomic Bose- Einstein 
condensate The measurements of the collective excitations have been found to be in excellent agreement with 
the theoretical predictions from mean- field theory for condensate fractions near unity ||- Such detailed studies 
of the collective excitations opens the door to investigating superfluid effects in atomic BECs, and in particular the 
general relation between BEC and superfluidity (see, for example, the article by Huang in Ref. (r|). 

In this paper we investigate the intereference of two BECs formed in small atomic samples composed of 10 3 — 10 6 
atoms, typical of current experiments. Our motivation for this theoretical investigation is that we previously showed 
that in small atomic samples the macroscopic wave function exhibits collapses and revivals, the first collapse occuring 
on a few seconds time scale jl4f| . In that work the notion of Bose broken symmetry |l5f| was invoked and the quantum 
state of the condensate was taken as a wave packet composed of states of fixed number of particles, so that, according 
to the uncertainty relationship between particle number and phase, the condensate phase can become well defined. 
The collapse of the macroscopic wave function arises from the fact that, due to many-body interactions, the chemical 
potential is different for each particle number present in the wave packet, and the revivals are a direct consequence 
of the discreteness of the particle number (lj]. Based on this theory one would naively expect that if two BECs 
are interfered the visibility of the interference pattern should also exhibit the collapses and revivals, hence providing 
a direct experimental signature of the effect. However, before jumping to conclusions it is important to note that 
the notion of Bose broken symmetry is only strictly applicable in the thermodynamic limit, where the system size 
and the number of particles tend to infinity with the density held constant, so that another approach should be 
utilized to verify the predicted collapse and revivals. Such an approach, freed from the thermodynamic limit and 
Bose broken symmetry arguments, has emerged in the last year p3|. Javanainen and Yoo |l6| first showed that 
an interference pattern between two condensates, and hence a relative phase between the two condensates, could 
arise from an explicit model of atomic detections: as atomic detections are performed an uncertainty in the number 
of atoms in each condensate is built up, since we don't know which condensate any atom is removed from. The 
uncertainty relationship between the relative atom number of the condendates and the relative phase then allows a 
well defined relative phase to emerge. 

The remainder of this paper is organized as follows: Section II gives the basic theory for BEC in a trapped gas of 
atoms as described by the Gross-Pitaevskii equation and within the Thomas-Fermi approximation. The interference 
between two condensates using the notion of Bose broken symmetry is described in Sec. Ill, and the collapse and 
revival times for the macroscopic wave function are evaluated for current experimental conditions. In Sec. IV we 
approach the interference of two condensates using a quantum anharmonic oscillator approximation to the full field- 
theoretical model of Sec. II, and the atomic detection scheme of Ref. fig] . Here we verify that the collapse and 
revivals are still present without explicitly invoking Bose broken symmetry. In Sec. V we briefly discuss the quantum 
anharmonic oscillator model for condensates described by coherent states, ie. assuming Bose broken symmetry, and 
demonstrate that the details of the collapses and revivals depends on the quantum state of the condensate. Finally, 
Sec. VI gives our summary and conclusions. 



II. BASIC THEORY 

In this Section we introduce the second quantized Hamiltonian for BEC in a trapped gas of atoms, and then discuss 
the Gross-Pitaevskii equation for the macroscopic wave function and the Thomas-Fermi approximation. Although 
these topics are amply covered elsewhere we include them for completeness and to clarify our notation. 

A. Second quantized Hamiltonian 

Our starting point is the second quantized Hamiltonian for a system of bosonic atoms confined in a trap potential 



H(t) = / dr 
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where ^(r, t) and 7p'{r,t) are the Heisenberg picture field operators which annihilate and create atoms at position r, 
and obey the equal-time commutation relation [V>(r, t),^'(r', t)] = S(r— r') appropriate to bosons, the (single-particle) 
trap potential is taken of the form [p5| p7| 

V(r) = ±m(u>° ± ) 2 (r 2 x + \ 2 z 2 ), (2) 

with r = (rj_, z), r± being the transverse position coordinate assuming cylindrical symmetry of the trap potential in 
the transverse plane, and z the longitudinal coordinate, m is the atomic mass, ujj_ the transverse angular frequency 
of the trap, A = lu®/lu^ is the ratio of the longitudinal to transverse frequencies, and Uq — ^Trfi 2 a/m measures the 
strength of the two-body interaction, a being the s-wave scattering length. Here we consider a repulsive interaction 
so that a > 0. 



B. Gross-Pitaevskii equation 



Here we consider the standard treatment of the condensate in a Bose gas |28| . At zero temperature, and for a 
weakly interacting gas, the particles may be assumed to be predominantly in the condensate. The assumption of zero 
temperature is not too restrictive since atomic condensates can be prepared with condensate fractions close to unity 
||[|. The Schrodinger equation for the state vector of the system is 



*ft^|*(f)>=F(0)|*(t)>, 



(3) 



and in the time-dependent Hartree approximation [E9J-B11 the state vector for a system of N particles is written as 



= (NW- 1 / 2 



dr^ w (r,t)^(r,0) 



JY 



|o>, 



(4) 



where ^>iv(r, t) is the normalized single-particle wave function, and |0) the vacuum state. Then proceeding in the 
usual manner |^,[5l| from the Schrodinger equation (||) we obtain the self-consistent nonlinear Schrodinger equation 
or Gross-Pitaevskii equation [E8[ generalized to include the magneto-optical harmonic trap p9 E7j 



ih 



dtfj 



N 



dt 
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N, 



(5) 



We are interested in the ground state solution of this equation for which we set ^jv(r, t) = exp(— i[iNt/h)(f>N(r), giving 
the stationary Gross-Pitaevskii equation 



HN<t> 



N 



^V 2 + imK) 2 (ri + A 2 z 2 ) + NU \M 2 
Am A 



(6) 



where the macroscopic wave function <fr]y(r) for the iV-particle system is normalized to unity. The ground-state energy 
of the system of TV atoms is derived from the Ginzburg-Pitaevskii-Gross energy functional |32| 



£ N = N I dr 
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and the parameter /in is given by /in = dS^/dN. By a slight extension of the usual terminology will be called the 
chemical potential of the N-particle condensate. The stationary Gross-Pitaevskii equation follows from the variational 
principle S(£n — N[in J dr\(f)N\ 2 ) = 0, and using this we find the further exact result 



, _ d/j, N 



U J dr|0jv| 



NU 



which can be rearranged to yield 



U Q J dr\<l> N \ 



N d 
1 + ~2dN 



(8) 



(9) 



which we shall use below. 



3 



C. Thomas-Fermi approximation 



Edwards and Burnett [ p5[ , and Ruprecht et. al. [£6| have solved Eq. (g) numerically to obtain both the condensate 
wave functions and the chemical potentials /in as functions of N. Here we follow Baym and Pethick p7j and use the 
Thomas-Fermi approximation for the macroscopic wave function, in which the kinetic energy term in the stationary 
Gross-Pitaevskii equation is neglected. This yields the approximate macroscopic wave function 



<Mr±, z) = (NU )- 1/2 ^n - imK) 2 (ri + A 2 z 2 ), (10) 

when the argument of the square root is greater than or equal to zero, zero otherwise. Requiring that the macroscopic 
wave function be normalized to unity yields the expression for the chemical potential 

where aj_ = (fi/mu') 1 '' 2 is the characteristic transverse length scale of the linear potential. The Thomas-Fermi 
approximation is valid if the coherence length r co h — (87ran)~ 1//2 p8[ , n being the mean density, is less than the 
characteristic size rjv of the atomic cloud. Following Baym and Pethick (27) we find that 

C =^=r^y 2/5 , as) 

r N \ a ± J 

which should be much less than unity for the Thomas-Fermi approximation to apply. In the first four columns of 
Table I we show the values of the parameters for some current experiments with a positive scattering length, and the 
fifth column shows the corresponding values of £, and we see that the Thomas-Fermi theory should be valid in all 
cases. 

The explicit form of the chemical potential in Eq. (11) can be used to evaluate the overlap integral in Eq. (||) 
giving 

,4 10 , 



U Q / dr|&v|=y/4, (13) 



a result we shall use later. 



III. QUANTUM INTERFERENCE WITH BOSE BROKEN SYMMETRY 

In this Section we investigate quantum interference between two BECs using the notion of Bose broken gauge 
symmetry p5[ . According to this notion the macroscopic wave function for each condensate attains a fixed but 
arbitrary phase as a result of the spontaneous breaking of U(l) gauge symmetry during the condensation process. 
The U(l) gauge symmetry is associated with particle number conservation, so that in the condensed system the 
particle number is not fixed. Then, as a result of the number-phase uncertainty relationship ANArj ss 1, the phase 
X] need not be indeterminate, and interference between two BECs is possible. Here we employ a description of the 
condensate state vector as a wavepacket of states of fixed particle number N in the condensate to apply the notion 
of Bose broken gauge symmetry. In contrast, in Section ^ we will show atomic detections establish a coherent 
wavepacket without the need to build in Bose broken symmetry from the outset. 



A. Wavepacket description 

The definition of the macroscopic wave function given in the last Section is rigorous in the thermodynamic limit 
p8fl . To investigate the case of a small condensate, say a few thousand atoms, we employ a wavepacket description. In 
particular, the wavepacket description is intended to reflect the quantum coherence of the condensate for measurement 



times short compared to relaxation times [22 33|]. In contrast, for measurements performed over times long compared 
to relaxation times the quantum coherence of the condensate is destroyed and the macroscopic wave function vanishes. 

Here we employ a wavepacket composed of states of fixed number of particles N in the condensate, with expansion 
coefficients = |a7v|e^ N , hence retaining quantum coherence. The present description of BEC, due to Barnett et. 



4 



al. j22[ , is therefore different from the conventional 77-ensemble which employs a wave-packet of states corresponding 
to different total number of particles, condensate plus non-condensate, and is generally only applicable in the ther- 
modynamic limit |33). A coherent state with ajv = N N / 2 e zNri e~ N / 2 / \/W. suggests itself, but the states associated 
with BEC do not generally possess such complete phase-coherence [Q. Nevertheless, a pure state description of 
the condensate may be rendered plausible as follows: Below the Einstein condensation temperature the many-body 
ground state of the system becomes macroscopically occupied, yielding a condensate which should be considered an 
open quantum system in contact via many-body interactions with the environment or reservoir composed of the non- 
condensate atoms. It is known from the work of Zurek et. al. |35) and Gallis that for a system, in their case a 
harmonic oscillator, in interaction with an environment certain pure states show considerable stability against loss of 
quantum coherence, and that in the weak coupling limit the pure states of maximal stability are the coherent states. 
Thus, the quantum coherence of the condensate may be reasonably represented by a pure state, though perhaps not 
precisely a coherent state since weak coupling may not hold. This argument does not depend on the size of the system, 
except that the non-condensate atoms may be viewed as a reservoir, and the conclusions therefore apply even for 
small condensates far removed from the thermodynamic limit. 

We further assume that the particle number distribution |ajv| 2 is sharply peaked with variance A A around a 
mean particle number A. As a concrete example we take a Poissonian distribution for which A A = A 1 / 2 , which is 
reasonable since the number distribution |ajv| 2 of the condensate should be approximately that of a coherent state, 
though the phases (jv may not be so correlated. Under these conditions we may expand the chemical potential /in 
around the mean number N in a Taylor series 

MJ v = Mjv + (N - AO/4 + \( N ~ N) 2 ^ + ■■■, (14) 

where 

If we compute the ratio of the third term to the second term in this expansion for A = N + AN, then we find that 
the magnitude of this ratio varies as N- 1 / 2 . Thus, even for a low particle number A = 10 3 the second term in the 
expansion is the dominant correction. In the remainder of this paper we shall therefore only retain the first two terms 
in the above expansion. 

B. Condensate state vector 

We are now equipped to construct the wavepacket for the condensate state vector, which we write as 

= 5^ a N e- lN ^ Nt/h {N\)- 1 ' 2 (tf) N |0), (16) 

AT 

where the exponential contains the term N/in since there are A particles each of chemical potential \x at, the operator 
a) = j d 3 r<j>N(r)ip^ (r, 0) creates particles with distribution </>jv(r) JUl) with [a,al] = 1, and |0) is the vacuum state. 
The approximations employed in Eq. (^) are tantamount to the Hartree approximation. In general the Schrodinger 
field annihilation operator can be written as a mode expansion over single-particle states as 

V>(r,Q) = ^2a a ip a (r) 

a 

= a^(r)+^(r,0), (17) 

where {^(r)} are a complete orthonormal basis set. Here in the last line we have chosen fo(r) = <^jv(r) as one 
member of the complete orthonormal set, and identified a — ciq. Then by construction the first term in the mode 
expansion acts only on the condensate state vector, whereas the second term ip(r,0) accounts for the non-condensate 
atoms. 



C. Collapses and revivals 

The state vector ( |l6| ) may be used to calculate quantum expectation values of the condensate. Of interest to us 
here are the one-particle reduced density matrix representing the condensate 
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p 1 (r,r' ) t) = (*(t)|StS|*(t))^(r)^(r') 

= ^(r)^(r'), (18) 



and the macroscopic wave function 



where 



(^(r,t)) = <*(t)|S|*(t))^(r) 

= N^^e-^^it), (19) 

^(t) = /f«w-i«A'[cos(2/ 1 ' J v(iV - N)t/K) - i sin(2/4(JV - N)t/h)}, (20) 

and n = [i^ + ^^'n * s ^ ne ne ^ cnem i ca l potential of the condensate. Using Eq. ( pi] ) for the chemical potential /zjv 
in the Thomas-Fermi approximation we obtain [i = |/ijy = £n/N, and the net chemical potential /i is equal to the 
mean energy per particle |27f| . In obtaining these results we have used the first two terms in the expansion of the 
chemical potential ([l4j). The one-particle reduced density matri x (|l8|) is of the classic factorized form representing 
the off-diagonal-long-range order (ODLRO) associated with BEC fl37| , [38| . This conclusion holds for any mean particle 
number, as long as the mean field approximations employed are valid. Here the ODLRO extends over separations 
|r — r'| w r^y, the spatial scale of the condensate. In the experiments the one-particle reduced density matrix is not 
measured but rather the momentum space (velocity) distribution is obtained by releasing the atoms from the trap, 
letting them fall under gravity, and imaging them. The momentum spread will then be AK 2ir/rjy, so that the 
ODLRO is transferred to a sharp spike in the imaged atomic distribution 

Turning now to the macroscopic wave function (119), we see that the factor Tjq{t) takes the form of a weighted sum 
of trigonometric functions with different frequencies. Such sums are well known from the Jaynes-Cummings model 
of quantum optics, which describes the interaction between a single-mode radiation field and a two-level atom, and 
give rise to the phenomenon of collapses and revivals |}9|, and the same is expected here. Collapse and revivals also 
appear in the relative phase between two superfluids or superconductors . Directly from the form of Eq. (^0|) we 
see that I-FjvWI ^ s periodic in time with period 

T R = nh/i/jf, (21) 

and the revivals occur with this period. The revivals result from the fact that the sum in Eq. (|2^) is over the discrete 
particle number, so that they are a direct result of the granularity of matter. The collapses depend on the choice of the 
number distribution |ajv| 2 . For our purposes we only need the variance A7V in particle number and the assumption 
that the phases are correlated enough that J-^(0) does not vanish exactly. Then J-ff(t) is periodic and has a maximum 
in magnitude at t = t max . At this time the net phases for each TV are such that they add most constructively in 
the sum in Eq. (pp[). As time increases these phase relations will initially be lost thus producing collapse of the 
magnitude of !Fjf(t) — > 0, until the system revives at t — t max + Tjv- The collapse time t co a may be estimated by 
looking at the spread of frequencies present in the wavepacket for particle numbers between N = N ± AN/2, which 
yields Ail = 2^,'^AN/h, and t co u sa 27r/Af2. Gathering our results together for the collapse and revival times we have 

The collapse phenomenon actually occurs under far more general conditions than reflected by the approximations 
used here, the essential ingredient being dispersion of the chemical potential /ijy over the particle number variance 
AN. Landau damping of plasma oscillations in an electron plasma is another example of decay of a coherent state. 
In contrast, the exactly periodic revivals arise from the linear dependence of /Ujv on (N — N) employed in Eq. ( |l4| 
). If higher-order corrections are retained in the expansion ( |i"4| ) the revivals are no longer perfectly periodic, and 
diminish with increasing time. We also note that according to Eq. (|l8|) the single-particle reduced density matrix is 
insensitive to the collapse and revivals of the macroscopic wave function, which is then also the case for the atomic 
BEC experiments |],§ . 



D. Thermodynamic limit 

The thermodynamic limit of these results must be taken with care. In particular, in order to maintain a constant 
density as the mean particle number is increased it is necessary to concomitantly increase the linear trap size as 
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r t cx TV 1 / 6 . Then cx N and t co u oc A 1 / 2 , so that collapse and revivals become irrelevant for N — > oo. In addition, 
we find that .Fjy(O) — ► & lr> by the following argument: The approximate uncertainty relation ANA-q w 1 holds for the 
number and phase fluctuations of the condensate. Then, as AA = A 1 / 2 — > oo we have A?] — > 0, which is the case 
for a coherent state with phases Cn = ^Vj?< Thus, in the thermodynamic limit the quantum state of the condensate 
approaches a coherent state for which we find ^jv(O) = e'" 7 . We then have {ip(r,t)) — e^A 1 / 2 ^ (r)e~ vt / ri , and this 
is precisely the limit in which the macroscopic wave function acts as an order parameter |28[ . 

E. Interference between two condensates 

We now consider the case that the atomic BEC is interfered with a second large condensate whose macroscopic 
wave function we write as C exp[i(k • r — fit/h + x)]- Here we have assumed that the second condensate is spatially 
large and write it as a plane-wave with wavevector k, and has the same chemical potential as the first for simplicity. 
We further assume that the second condensate is composed of a large number of atoms so that collapses and revivals 
are not an issue, the constant C is real, and \ is t ne arbitrary but fixed phase of the condensate. Then, when the 
two condensates are made to interfere, say by dropping them on top of each other, the resulting interference pattern 
is described by the cross term 

T(r,i) = C(ijj(r,t)e- l{k - r+x) + ^(r, i)e l(k ' r+x) ) 

= 2CN 1 / 2 ^(r)\^(t)\ cos(k • r + X - arg(T N {t))). (23) 

In the thermodynamic limit where Ffj{t) — > e %T> this interference pattern takes the form of a stationary spatial 
modulation of wavevector k over the spatial extent of the condensate (we assume |k|r^ >> 1 so that an interference 
pattern is evident). However, due to the factor appearing in the interference pattern it is clear that the 

interference will display the collapses and revivals predicted for the macroscopic wave function: when the macroscopic 
wave function collapses and revives so does the visibility of the interference pattern with the same period Tpj. 

Thus quantum interference of two condensates provides an experimental signature of the predicted collapses and 
revivals. Furthermore, the arguments are easily extended to the case that both condensates display collapses and 
revivals, that is two identical condensates, and the same conclusions apply. 

F. Experimental parameters 

The collapse and revival times calculated for the current experiments are shown in columns six and seven of Table 
I, respectively. Here we see that except for the weak trap in Rb (second row) the collapse times are all less than one 
second, indicating that collapse of the macroscopic wave function occurs in these traps. In contrast, the revival times 
are considerably longer. However, the 6 s revival time found for the strong trap in Rb is still within the condensate 
lifetime of 15 s quoted for that experiment jj], though the effects of dissipation may quench the revival. In this 
respect, we note that the recent Rb experiment (third row) with a higher trap frequency and mean particle number 
has a lower revival time of 3.75 s. 



IV. QUANTUM PHASE FROM MEASUREMENTS 

In this Section we look at the build up of quantum coherence between two Bose-Einstein condensates initially in 
Fock states. Since we start from Fock states no phase is initially assumed with a phase only established via the 
measurement process. The system we consider was first proposed by Javanainen and Yoo fl6|| , and consists of two 
Bose-Einstein condensates which are dropped on top of one another. Each of the condensates consist of n atoms with 
momentum k\ and &2 directed along the x-axis, respectively. Atoms are detected on a screen placed below the two 
condensates. A detection at some position x is represented by the field operator for the sum of the two condensates 

*(x) = (ai + &2e*+to) , (24) 

where <f>{x) = {k% — ki)x, and hi and a>2 are the atom annihilation operators for the first and second Bose-Einstein 
condensates respectively (see below). An interference pattern is generated from the two condensates because every 
detection of an atom introduces uncertainty into the atom number in each condensate since we do not know which 
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condensate this atom came from. However, the total number in both condensates is always known since it is just 
the initial total number minus the number of detections we have observed. The uncertainty relationship between the 
relative atom number and phase between these condensates allows us to establish a relative phase between them to 
some precision. 

The effect of collisions on the establishment of quantum phase has been studied previously [^lj. Here we extend 
this work to demonstrate collapses and revivals in the visibility pattern between two interfering condensates. 



A. Single-mode approximation 



For the purposes of developing a model based on atomic detection the full quantum field theory involving the field 
operators is cumbersome. Thus, here we reduce the full problem by using a mode expansion and truncating. In 
particular, we write the Heisenberg field annihilation operator as a mode expansion over single-particle states as 



-ifi a t/h 



= a(t)<f >R (T)e- i »* t / h + i,(r,t), 



(25) 



where {(p a (r)} are a complete orthonormal basis set, and {fi a } the corresponding eigenvalues. Here in the last line we 
have chosen tpo( r ) = ^jvl 1 *) as one member of the complete orthonormal set, and identified /xq = A* at and a — do- Then 
by construction the first term in the mode expansion acts only on the condensate state vector, whereas the second 
term i/j(r,t) accounts for the non-condensate atoms. Then substituting the mode expansion in the second-quantized 
Hamiltonian (^), retaining only the first term representing the condensate, and using the Gross-Pitaevskii Eq. (|^), 
we obtain the following single-mode quantum Hamiltonian for the condensate in the Schrodinger picture 



H(0) 



a' ae 



N 



a)a)aa, 



(26) 



where 



fin = Uq / dr\(f> 



,J N\ 



(27) 



which using Eq. ( |13| ) yields fin = -^A 4 ^, K being the collision rate between condensate atoms. The first term on the 
right-hand-side, which gives rise to an energy shift, is the number operator a' a times the single-particle contribution 
to the energy per particle, 



-N 



dr 



2m 



iv^p + nr)!^! 2 



(28) 



whereas the second term accounts for many-body interactions. To proceed we introduce an interaction picture defined 
by the transformation 



Hj(t) = UHt) 



H s - ih 



dU(t) 
dt 



U(t), U(t) 



^—i(efj-~hK 1 /2)a J *at/h 



which yields the single-mode Hamiltonian 



H{t) = ^a)\ 



(29) 



(30) 



where we have dropped the subscript I on the interaction picture operators for simplicity in notation. In the following 
Sections we shall use this single-mode Hamiltonian, which is that of a quantum anharmonic oscillator, as a model for 
an individual atomic BEC. 

A basic assumption underlying this model is that the quantum state describing the condensate involves only particle 
numbers n ~ N, that is the particle number remains close to the mean number. The single-mode Hamiltonian ( |30| ) 
gives rise to a particle number dependent phase-shift \n = K N = ^{ifqN/Ti (see Sec. V). However, Eq. ( pc| ) reveals 
a particle number dependent phase-shift xn = 2^j,^N/Tl. The discrepancy between these two results arises since 
the phase-shift in the single-mode theory is an average over the spatial profile of the macroscopic wave function, as 
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evidenced by the overlap integral in the definition of the collision rate k in Eq. (g7|) . This averaged value is less than 
the peak value of the actual phase-shift which appears in expression ((2^). In the following we shall replace the factor 
t£ by 2 in the definition of k 

Tin = 2fjhff, (31) 

so that the correct phase-shift is produced by the quantum anharmonic oscillator model. This relation provides 
the connection between the single-mode theory parameters and those obtained from the treatment based on the 
Gross-Pitaevskii equation. 

In general we shall consider two interfering condensates. For the i th condensate described by single-mode operators 
a,i,a\ we assign a wave vector ki, directed along the x-axis for simplicity, so that the interference between the 
condensates can be discussed. The resulting interference pattern would be further modulated by the detailed spatial 
structure of each condensate, as determined by the solution of the Gross-Pitaevskii equation, but this is of secondary 
importance and we do not discuss it here. In general, we only assume that any interference is well resolved within the 
spatial extent of the condensates. 



B. Effect of Collisions 



The effect of collisions can be numerically simulated via a Monte Carlo wave function method. The Monte Carlo 
method is used to simulate both the detection process of rate 7 (as quantum jumps) and the time evolutions between 
detections (as evolution of the wave function between the jumps). The initial Fock state, say \n,n), becomes an 
expansion of to + 1 entangled Fock states after m detections. Let us write the state vector after to detections as 



\<p m ) = X! Cfe l n - m + k,n - k), 



(32) 



k=0 



which is normalized so that X) l c fc| = 1. The effects of the collisions are included in the time evolution operator for 
the two condensates 



U (t) = e - i ( ffl + fra )*/ R = exp 1 - — 



a\a 2 



(33) 



with Hi the Hamiltonian for the i th condensate. Thus after m detections we have prepared a state \<p m ) which then 
experiences the above evolution operator for a time t followed by a further detection. The conditional visibility of 
this detection is 



p (x|xi , . . . , x nl: 
The expectation value for the total number operator is 

(ip m \U^ (t) (a{ai + a\a^j U (t) \(p m ) = 2n - to, 



(34) 



(35) 



and, as expected, it is just the total initial number of atoms minus the number of detections. The expectation value 
of the cross term which gives the size of the conditional visibility is 



(<p m \tf [t)a\a 2 U (t)\^ m )e^ 

m 

= c* k Ck-i\/ (n — k + 1) (n — to + k) exp{ihd [2k — (to + 1) + i<j> (x)]} . 



k=l 



Putting Eqs. (||) and (||) into Eq. (||) we obtain 

m 

p (x\xi, . . . , x m , t) oc n — m/2 + A (k) cos [(f) (x) + nt (2k — m — 1) — 6&] , 



(36) 



(37) 



k=l 



where we have defined the phase 0& from the coefficients of the state vector as c^Ck-i = A^e k and the weighting 
function A is 
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A (k) = A k ^(n- k + l)(n-m + k). 



(38) 



We take Ofc outside the summation since it is defined as the relative phase between neighboring number states which 
we will assume to be fairly constant for large m. As we detect more and more atoms the resulting entangled state 
approaches something that resembles a coherent state where the relative phase between neighboring number states 
(in a number state basis) is identical to the phase of this state. Let this fairly constant relative phase be denoted 
as O which is a good estimate of the relative phase between the two condensates. By expanding the cosine terms in 
Eq. ( |37| ) to separate the time dependence from the phase terms, and noting that the resulting summation over sine 
functions vanishes due to the cancelation of positive and negative frequency components, we obtain 

m 

p (x\x%, . . . , x m , t) oc n — m/2 + A (k) cos [(2k — m — 1) nt] cos [<f> (x) — 0] . (39) 

fe=i 

The visibility of the interference pattern is therefore determined from a weighted sum over cosines of differing fre- 
quencies. These frequencies depends on the parity of the expression enclosed by the round brackets, 2k — m — 1. 
When this is odd (m is even) the frequencies are {1 — m, . . . , —3, — 1, 1, 3, . . . , m — 1} whereas when it is even (m is 
odd) they are {1 — m, . . . , —2, 0, 2, . . . , m — 1} . In both cases the revival period of the visibility is 

T=- = —, (40) 

with N — n. This revival period is precisely half that obtained assuming Bose broken symmetry in Sec. Ill, and this 
difference shall be taken up in Sec. V. For the case where m is even so that we have a sum over odd frequencies, 
the cosine term in Eq. ( J3S| ) alternates between plus or minus 1 at subsequent revival times. Since the visibility is by 
definition a positive quantity, this alternate sign change in the cosine term represents a ir phase shift of the interference 
pattern at alternate revivals. 

For collision to detection rate ratio of one we expect from Eq. ( [40| ) a period of n when time is measured in units 
of the inverse detection rate 7. This agrees very well to the collapse and revival period displayed in the middle graph 
of Fig. (0). This figure and the ir phase shift will be explained in detail in Sec. V. 



C. Numerical results 



In a numerical simulation the effects of collisions between subsequent detections can be easily modeled with the use 
of the Monte Carlo wave function method §42] . We use the following effective Hamiltonian 



H eff - TT 



(a\a-, 



(41) 



where the detection and collision rates are 7 and k respectively. The detections are then turned off and the system 
undergoes coherent evolution by the Hamiltonian 



2 



(42) 



We have modeled the collisions that occur only between the atoms of the same condensate. A cross-collision term 
(a\a\a\a2) between the condensates is not included. The size of the coefficient of this term depends on the physical 
geometry of the situation (the overlap of the two condensates) ranging from zero to hn. Setting this coefficient to 
zero we are taking the worst case scenario where the effects of the collisions are the strongest. Alternatively, setting 
the cross-collison term to fik completes the square in the Hamiltonian so that the subsequent evolution depends only 
on the total atom number. The effect of the collisions would then be to rotate the phase of the entire state vector, 
and the coherence between the individual entangled number states would be preserved. 

In each run of the numerical simulation the state vector experiences three different regimes. Initially a sequence 
of detections are accumulated to prepare the entangled state, after which the detections are turned off. During the 
coherent evolution stage free of detections, the conditional visibility undergoes collapses and revivals due to collisions. 
Finally, the detections are turned on again. If the detections are turned on when the visibility is zero, during a 
collapse, it is quickly re-established by the second sequence of detections. If, however, the detections are turn on 
when the visibility is in a revival phase the visibility starts from this non-zero value and then quickly increases to 
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one. This behavior is seen in Fig.|T] where the detection process occurs at a much faster rate than the collapse and 
revivals, and we thus graph the single run on three separate set of axes: The top graph shows the initial growth of 
the conditional visibility due to the detections of atoms from the condensates, and the visibility quickly increases 
to a value close to unity after 100 atoms arc detected. (In these simulations the total number of atoms is 10 5 and 
200 atomic detections are made, so the assumption underlying the quantum anharmonic oscillator model that the 
mean number of atoms varies little is well obeyed.) The middle graph shows what happens to the visibility once 
the detection process is turned off and the state vector undergoes coherent evolution due to the Hamiltonian in Eq. 
(ff2|), and collapses and revivals of the visibility are clearly evident. This behavior is reminiscent of the collapses and 
revivals in the Jaynes-Cummings model of quantum optics for a two-photon process which also displays the periodic 
revivals of Fig.[l]. The bottom graph displays the subsequent evolution of the visibility when the detection process 
is turned on again. Due to the collapses and revivals during the coherent evolution, the initial visibility for the final 
stage depends on the time at which the detections are reinitiated. Whatever the initial visibility, this second detection 
proceeds very rapidly to increase the visibility to a value close to unity. 

The cleanliness and exhibition of full revivals in the middle graph of Fig.l can be understood by noticing that 
the two condensate system is a closed system (no loss mechanism) undergoing coherent evolution. The state vector 
after the initial sequence of detections is an expansion of entangled number states. The effect of collisions during the 
coherent evolution is to rotate the phase of the coefficients of each entangled state by an amount proportional to the 
sum of the squares of the number of atoms in each condensate. Thus the phase of the coefficients of this state vector 
rotate at differing frequencies. 

So far we have shown collapses and revivals in the conditional visibility of the interference pattern. For something 
more relevant to an experimental situation we would like to look at variables associated with the actual observed 
interference patterns. The phase shift of an interference pattern is a direct measure of the relative phase established 
between the two condensates. Thus, let us consider the following procedure: Firstly, we prepare a state vector of the 
two condensates with an established relative phase via measurements, and consider this entangled state between the 
two condensates after the detections to be our initial state which possesses some degree of coherence. This state is 
then allowed to undergo coherent evolution with no detections, and finally we turn on the measurement process after 
an elapsed time and collect our second sequence of measurements. The phase of the resulting interference pattern 
is calculated with respect to the phase of the interference pattern we observed previously from the first sequence of 
measurements. Now we re-prepare the initial state and repeat the second sequence of measurements and subsequent 
calculation of phase. Repeating this many times we obtain a set of relative phases between the first and second set 
of measurements, the idea being that if the time elapsed between the two detection regimes corresponds to some 
multiple of the full revival time then this set of relative phases should be sharply peaked at zero. For other elapsed 
times we may expect to see partial revivals. We do not need to numerically calculate the coherent evolution but what 
we need from the numerical simulation is the coefficients of the prepared state \(p m )- Its coherent evolution due to the 
Hamiltonian H describing the collisions previously given by Eq. ( [42| ) is 

\<Pm (*)) = ex P (-iHt/k) \ip m ) 
i 



ex P I 



n 2 + (n — mf 

x Cfc exp (— in [— mk + fc 2 ] t) \n — m + k, n — k) . (43) 

fe=0 

We use the phase eigenstates for the atom number difference between the two condensates 



£ exp 

n 1 ,n 2 



- (m - n 2 )(f> 



K,n 2 ), (44) 



which has a factor of one- half in the exponential since the state \ipm) has a fixed total atom number of 2n — m so that 
the atom number difference [n\ — n%) is quantized in units of 2. Thus this factor is required so that <f> is the relative 
phase between the condensates. The probability distribution of the phase 4> after elapsed time t is 



Cfc exp (— i [k (k — m) nt + 2k(j>\) 



k=0 



(45) 



This probability is a function of two variables, phase and elapsed time, and is evaluated numerically. We display the 
probability distributions in "birds-eye" plots via the "image" command using MATLAB© . Figure (Ha) displays the 



11 



probability distribution when we have made an even number of detections (m even) . The white regions denotes the 
peaks with the black ones corresponding to the valleys. We see a sharp peak at zero time about zero phase difference 
between the two interference patterns, this is not surprising since no coherent evolution has occurred. By the time 
we have evolved for one revival period we obtain a sharp peak at n radians corresponding to the first revival time 
with 7r phase shift due to detecting an even number of atoms described in the previous section. The second revival 
occurs at a phase difference of zero radians as predicted. Away from these revival times, the phase distribution is 
not flat but displays many partial revivals. Figure (Qb) is a zoomed in view of Fig. (||a) between 0.2 and 0.5 revival 
periods, we can clearly see the partial revival at 0.2 which consists of five peaks. In fact we see partial revivals at 
every integer fraction of a revival period provided the resolution is good enough. We illustrate this by labelling to 
the left of the graph with the appropriate integer fractions corresponding to the particular partial revival. Figure (||) 
shows the distribution for the other case where the number of detections is odd. As predicted there is no it phase 
shifts. Again, we see partial revivals when we zoom in between 0.2 and 0.5 revival periods as shown in Fig. (^b). 



V. QUANTUM OSCILLATOR MODEL WITH BROKEN SYMMETRY 



We will show in this brief Section that collapses and revivals also arise for the quantum anharmonic oscillator model 
for initial coherent states instead of preparing an entangled state from an initial detection process. In the language of 
the previous Section we are considering the interference between two coherent states including the effects of collisions. 
Thus we treat the Bose-Einstein condensates as coherent states, analogous to the treatment in Sec. Ill, imposing on 
them a relative phase whereas previously we establish this phase via measurements. 

The Hamiltonian for the two condensates is 



\ ^ 2 

H = - ^ M (otoi) , 



(46) 



giving the Heisenberg equation of motion for each field anihillation operator 

By inspection the time dependence of the on operator is 

hi (t) = cxp — — ^1 + 2a|aiJ Kt 

which yields the Heisenberg picture field operator for the sum of the two modes 

1 



- + 2a\ai ) Kt 



CL\ + CXp 



^1 + 2a\a,2*j Kt 



«2 



(47) 



(48) 



(49) 



where we have suppressed the spatial dependence. This yields the operator for the intensity of the atomic pattern as 

1 



& (t) & (t) = — |aj&i + a|^2 + a\ exp i (a\ai — a\a,2*j Kt a,2 + h.c.j 



(50) 



Thus, if we define the initial coherent states as \a) and \/3) for the modes Siand a-i respectively, the intensity is 
evaluated to be 



lex (a,/3\& (t)$(t)\a,/3) 



a*/3exp 



1)1' 



+(3*a exp 

with the aid of the following identity p3| 



(e~ iK 4 - 1) H 2 + (e iKt - 1) |/3| 2 ] } , 



(51) 



a e xa a a) = exp (e - * - 1) |a| 
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In the case of maximum visibility of the interference pattern we need the modes to be equal in amplitude, thus we 
set (3 = ae~ l( ^ where <fi is the phase difference between the modes. Substituting this into Eq. (51) we obtain 



lex |a| 2 { 



1 + exp 



a 



(COS Kt — 1) 



(52) 



The characteristics of collapses and revivals can be clearly seen in Eq. (^2|). Revivals occur at times t — 2im/ n (where 
n is an integer) with the shape of the collapses described by the exponential term whenever t is no longer a multiple 
of 2tt/k. The period of these revivals is 



27r nh 



(53) 



which is identical to Eq. ( pi|) with N — \a\ 2 . This gives twice the period of the previous case of revivals from a 
detection process (see Eq. (jfO|)), even though we have used identical Hamiltonians. This difference arises because we 
have two independent coherent states here whereas in the case of detecting atoms from two initial Fock states the 
total atom number is always fixed at the initial total number minus the number of atomic detections. This additional 
constraint on the total number means that the difference atom number operator inside the exponential in Eq. ( |50| ) 
is quantized in units of 2. The above expression for the period is only applicable when we are free to superpose total 
atom numbers, as in the case of two coherent states for which the difference atom number operator is quantized 
in units of 1. This gives a factor of 2 difference in the revival times between assuming an initial relative phase and 
establishing this phase in the dynamics of the revivals. The visibility, as described by Eq. ((5^), smoothly drops to 
its minimal value halfway between subsequent revivals, exactly where the state induced by detection would have an 
additional revival. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we have shown that within the approximation of Bose broken symmetry the macroscopic wave function 
in small atomic samples exhibits collapse and revivals in time while the BEC is maintained in the form of ODLRO. 
For current experiments the collapse time is a second or less. The revival time is longer but our results show that it 
may still fall within the condensate lifetime for some experiments. To detect the collapses and revivals experimentally 
a scheme is required which is sensitive to the macroscopic wave function directly, and this does not seem to be 
the case for the coherent light scattering methods previously discussed Q- However, Imamoglu and Kennedy 
p7[ and Javanainen |48| ] have recently proposed light scattering schemes involving two independent condensates 
coupled by a common excited state. These schemes rely on the fact that when one condensate is driven optically 
the light scattered from the other condensate has a non-zero value of the electric field and a phase proportional to 
the relative phase of the two condensates. By driving both wells and adjusting the phase difference of the fields 
the scattering can be suppressed via quantum interference, and this in turn determines the phase difference between 
the two condensates. The scattering is therefore sensitive to Bose broken gauge symmetry In addition, the light 
scattering rate is proportional to the magnitude of the macroscopic wave function, so these schemes could be used to 
detect the collapses and revivals experimentally. Alternatively, as we have shown here, a direct measurement of the 
collapse and revivals in the visibilty of the interference between two condensates can be made. 

In a second approach we have also studied the establishment of a relative phase between two interfering condensates 
using the explicit measurement model first proposed by Javanainen and Yoo p6| . This model is free from any 
assumptions concerning Bose broken symmetry or the thermodynamic limit, and is therefore applicable to small 
atomic samples. In this case we consider the two interfering condensates to be in number states initially, an extreme 
example for which there is no relative phase before atomic detections. These condensates are then prepared into 
an entangled state vector composed of entangled number eigenstates via the measurement process. The effects of 
collisions in the time evolution is to rotate the phases of the individual entangled eigenstates of the state vector. We 
observe the collapse and revivals of this state under coherent evolution when the measurement process is turned off. 
Accurate predictions of the period of these collapses and revivals were obtained. The simple anharmonic model of 
interference between two condensates also displays collapses and revivals. The period of these revivals is twice the 
time required for the previous case since we have neglected the constraint of the total atom number being fixed. 



ACKNOWLEDGEMENTS 



Dan Walls acknowledges support from the Office of Naval Research, the New Zealand Foundation for Research 
Science and Technology, and the Marsden Fund Useful discussions with S. Barnett, K. Burnett, and J. C. Garrison 



13 



are appreciated. 



14 



Expt. 




A 


a 


N 


c 


tcoll 


Tf) 


1 


471 rad s" 1 


2.8 


5.2 nm 


2000 


0.1 


0.13 s 


6 s 


2 


51 rad s _i 


2.8 


5.2 nm 


2000 


0.15 


2 s 


86 s 


3 


828 rad s" 1 


2.8 


5.2 nm 


4500 


0.05 


55 ms 


3.75 s 


4 


2010 rad s _i 


0.056 


4.9 nm 


5 x 10 b 


0.02 


0.32 s 


710 s 



1 . " Strong trap" experiment in Rb from Ref. 

2. "Weak trap" experiment in Rb from Ref. 

3. Experiment in Rb from Ref. 

4. Experiment in Na from Ref. 



I 



15 



M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Science 269, 198 (1995). 

K. B. Davis, M. O. Mewes, M. R. Andrews, N. J. van Druten, D. S. Durfee, D. M. Kurn, and W. Ketterle, Phys. Rev. 

Lett. 75, 3969 (1995). 

M. O. Mewes, M. R. Andrews, N. J. van Druten, D. M. Kurn, D. S. Durfee, and W. Ketterle, Phys. Rev. Lett. 77, 416 
(1996). 

M. O. Mewes, M. R. Andrews, N. J. van Druten, D. S. Durfee, D. M. Kurn, and W. Ketterle, Science 273, 84 (1996). 
D. S. Jin, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Phys. Rev. Lett. 77, 420 (1996). 
M. O. Mewes, M. R. Andrews, N. J. van Druten, D. M. Kurn, D. S. Durfee, C. G. Townsend, and W. Ketterle, Phys. Rev. 
Lett. 77, 988 (1996). 

D. S. Jin, M. R. Matthews, J. R. Ensher, C. E. Wieman, and E. A. Cornell, "Temperature-dependent damping and 
frequency shifts in collective excitations of a dilute Bose-Einstein condensate," submitted to Phys. Rev. Lett. (1996). 

M. O. Mewes, M. R. Andrews, D. M. Kurn, C. G. Townsend, and W. Ketterle, "An output coupler for Bose condensed 
atoms," submitted to Phys. Rev. Lett. (1996). 
A. L. Fetter, Phys. Rev. A 53, 4245 (1996). 

K. G. Singh and D. S. Rokshar, Phys. Rev. Lett. 77, 2360 (1996). 

M. Edwards, P. A. Ruprecht, K. Burnett, R. J. Dodd, and C. W. Clark, Phys. Rev. Lett. 77, 1671 (1996). 
S. Stringari, Phys. Rev. Lett. 77, 2360 (1996). 

A. Griffin, D. W. Snoke, and S. Stringari, Eds., Bose-Einstein Condensation (Cambridge University Press, Cambridge, 
1995). 

E. M. Wright, D. F. Walls, and J. C. Garrison, Phys. Rev. Lett. 77, 2158 (1996). 

For a comprehensive discussion of Bose broken symmetry see, for example, A. Griffin, Excitations in a Bose- condensed 
liquid (Cambridge University Press, Cambridge, 1993) Chap. 3. 
J. Javanainen and S. M. Yoo, Phys. Rev. Lett., 76 ,161 (1996). 

M. Naraschewski, H. Wallis, A. Schenzle, J. I. Cirac and P. Zoller, Phys. Rev. A 54 2185 (1996). 
Y. Castin and J. Dalibard (unpublished). 

J. I. Cirac, C. W. Gardiner, M. Naraschewski, and P. Zoller, Phys. Rev. A 54, 3714R (1996). 
T. Wong, M. J. Collett, and D. F. Walls, Phys. Rev. A 54 3718R (1996). 
T. Wong, M. J. Collett, S. M. Tan, and D. F. Walls (unpublished). 

S. M. Barnett, K. Burnett, and J. Vaccarro, "Why a condensate can be thought of as having a definite phase," to appear 
in the Feature Issue on Bose-Einstein Condensation in the Journal of Research of the National Institute of Standards and 
Technology (1996). 
K. Maimer, (preprint). 

K. Huang, Statistical mechanics (John Wiley, New York, 1987, 2nd Ed.). 
M. Edwards and K. Burnett, Phys. Rev. A 51, 1382 (1995). 

P. A. Ruprecht, M. J. Holland, K. Burnett, and M. Edwards, Phys. Rev. A 51, 4704 (1995). 
G. Baym and C. J. Pethick, Phys. Rev. Lett. 76, 6 (1996). 

E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, Part 2 (Pergamon Press, Oxford, 1989)pp. 85-118. 
P. A. M. Dirac, Proc. Camb. Philos. Soc. 26, 376 (1930). 

A. K. Kerman and S. E. Koonin, Annals of Physics 100, 332 (1976). 

B. Yoon and J. W. Negele, Phys. Rev. A 16, 1451 (1977). 

V. L. Ginzburg and L. P. Pitaevskii, Zh. Eksp. Teor. Fiz. 34, 1240 (1958) [Sov. Phys. JETP 7, 858 (1958)]; E. P. Gross, 
J. Math. Phys. 4, 195 (1963). 

P. C. Hohenberg and P. C. Martin, Ann. Phys. (New York) 34, 291 (1965). 

P. W. Anderson, Basic notions of condensed matter physics (Addison- Wesley, Reading, 1984)pp. 229-261. 

W. H. Zurek, S. Habib, and J. Pablo Paz, Phys. Rev. Lett. 70, 1187 (1993). 

M. R. Gallis, Phys. Rev. A 53, 655 (1996). 

O. Penrose and L. Onsager, Phys. Rev. 104, 576 (1956) 

C. N. Yang, Rev. Mod. Phys. 34, 694 (1962). 

The theoretical description of collapse and revivals in the Jaynes-Cummings model is described in J. H. Eberly, N. B. 
Narozhny, and J. J. Sanchez-Mondragon, Phys. Rev. Lett. 44, 1323 (1980). 

F. Sols, Physica B 194-196, 1389 (1994). 

T. Wong, M.J. Collett and D.F. Walls, to be published in Phys. Rev. A. 

For details on effective Hamiltonians in Monte Carlo methods see CW. Gardiner, A.S. Parkins and P. Zoller, Phys. Rev. A 
46, 4363 (1992); J. Dalibard, Y. Castin and K. Maimer, J. Opt. Soc. Am 10, (1993); H.J. Carmichael, An Open Systems 
Approach to Quantum Optics, Lecture Notes in Physics ml8 (Springer- Verlag, Berlin, 1993). 



1G 



[43] W.H. Louisell, Quantum Statistical Properties of Radiation, (Wiley, New York 1973) pp. 156. 
[44] J. Javanainen, Phys. Rev. Lett. 72, 2375 (1994). 

[45] L. You, M. Lewenstein, and J. Cooper, Phys. Rev. A 50, R3565 (1994). 
[46] J. Javanainen and J. Ruostekoski, Phys. Rev. A 52, 3033 (1995). 

[47] A. Imamoglu and T. A. B. Kennedy, "Optical signatures of Bose-broken symmetry," submitted to Phys. Rev. A (1996). 
[48] J. Javanainen (submitted to Phys. Rev. Lett.). 

FIG. 1. The top graph shows the growth in the conditional visibility as a function of time corresponding to 100 atomic 
detections. The middle displays the collapse and revivals of this conditional visibility when the detection process is turned off 
and the bottom graph shows the growth of this visibility when the detection process is turned on again. The total number of 
atoms in the two condensates was 10, 000 and we have used a collision to detection rate ratio of one (k = 7) for the initial and 
final detection periods. 

FIG. 2. Plot of the phase distribution when the number of detections (m) is even as a function of the turning on time. The 
phase is the relative phase between the first and second sequence of measurements in units of radians while the turning on time 
is in units of the revival period. The brightness of a region corresponds to the relative probability of obtaining a particular 
phase for a particular turning on time. The bright regions denotes peaks while the darker ones correspond to valleys. The 
entire plot is shown in (a) with a zoomed in plot between 0.2 to 0.5 times displayed in (b). 

FIG. 3. Plot of the phase distribution when the number of detections (m) is odd as a function of the turning on time. The 
entire plot is shown in (a) with a zoomed in plot between 0.2 to 0.5 times displayed in (b). 
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